Evidence for the Onset of Deconfinement 
from Longitudinal Momentum Distributions? 
Observation of the Softest Point of the Equation of State 
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We analyze longitudinal pion spectra from i?i a b = 2A GeV to ^/snn = 200 GeV within Landau's 
hydrodynamical model. From the measured data on the widths of the pion rapidity spectra, we 
extract the sound velocity c 2 s in the early stage of the reactions. It is found that the sound velocity 
has a local minimum (indicating a softest point in the equation of state, EoS) at -Ebeam = 30A GeV. 
This softening of the EoS is compatible with the assumption of the formation of a mixed phase at 
the onset of deconfinement. 



Over the last years, a wealth of detailed data in the 
20A — 160A GeV energy regime has become available. 
The systematic study of these data revealed surpris- 
ing (non-monotonous) structures in various observables 
around 30A GeV beam energy. Most notable irregular 
structures in that energy regime include, 

• the sharp maximum in the K+/7T+ ratio 0,0, 

• a step in the transverse momentum excitation func- 
tion (as seen through (mi) — too ) 00 

• an apparent change in the pion per participant ratio 
@ and 

• increased ratio fluctuations (due to missing data at 
low energies it is unknown if this is a local max- 
imum or an ongoing increase of the fluctuations) 



It has been speculated, that these observation hint to- 
wards the onset of deconfinement already at 30A GeV 
beam energy. Indeed, increased strangeness production 
Q and enhanced fluctuations have long been predicted 
as a sign of QGP formation H 0, S H Hfl El within dif- 
ferent frameworks and observables. The suggestion of an 
enhanced strangeness to entropy ratio (~ K/tt) as indica- 
tor for the onset of QGP formation was especially advo- 
cated in [l^. Also the high and approximately constant 
K ± inverse slopes of the tot spectra above ~ 30A GeV 
- the 'step' - was also found to be consistent with the as- 
sumption of a parton <-> hadron phase transition at low 
SPS energies 0,E1- Surprisingly, transport simulations 
(supplemented by recent lattice QCD (1QCD) calcula- 
tions) have also suggested that partonic degrees of free- 
dom mi ght already lead to visible effects at ~ 30A GeV 
[TH [psL Tl7|. Finally, the comparison of the thermody- 
namic parameters T and fin extracted from the transport 
models in the central overlap region [Tsj with the experi- 
mental systematics on chemical freeze-out configurations 
[I9II20II21I I in the T—/j,b plane do also suggest that a first 
glimpse on a deconfined state might be possible around 
10A - 30A GeV. 

In this letter, we explore whether similar irregu- 
larities are also present in the excitation function of 
longitudinal observables, namely rapidity distributions. 



Here we will employ Landau's hydrodynamical model 
IH El El 111 IH H3. This model entered the 
focus again after the most remarkable observation that 
the rapidity distributions at all investigated energies can 
be well described by a single Gaussian at each energy. 
The energy dependence of the width can also be rea- 
sonably described by the same model. For recent appli- 
cations of Landau's model to relativistic hadron-hadron 
and nucleus- nucleus interactions the reader is referred to 

mill El Him ( andRefs - therein). 

The main physics assumptions of Landau's picture are 
as follows: The collision of two Lorentz-contracted nuclei 
leads to full thermalization in a volume of size V/ ^/snn- 
This justifies the use of thermodynamics and establishes 
the system size and energy dependence. Usually, a simple 
equation of state p = c 2 s e with c 2 s = 1/3 (c s denotes the 
speed of sound) is assumed. For simplicity, chemical po- 
tentials are not taken into account. From these assump- 
tions follows a universal formula for the distribution of 
the produced entropy, determined mainly by the initial 
Lorentz contraction and Gaussian rapidity spectrum for 
newly produced particles. Under the condition that c s is 
independent of temperature, the rapidity density is given 
byHilH: 
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where if is a normalisation factor and m p is the proton 
mass. The model relates the observed particle multiplic- 
ity and distribution in a simple and direct way to the 
parameters of the QCD matter under consideration. 

Let us now analyze the available experimental data 
on rapidity distributions of negatively charged pions in 
terms of the Landau model. Fig. ^ shows the measured 
root mean square a y of the rapidity distribution of neg- 
atively charged pions in central Pb+Pb (Au+Au) reac- 
tions as a function of the beam rapidity. The dotted line 
indicates the Landau model predictions with the com- 
monly used constant sound velocity c 2 s — 1/3. The full 
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FIG. 1: The root mean square a y of the rapidity distributions 
of negatively charged pions in central Pb+Pb (Au+Au) reac- 
tions as a function of the beam rapidity y p . The dotted line 
indicates the Landau model prediction with = 1/3, while 
the full line shows a linear fit thr oug h the data points. Data 
(full symbols) are taken from F3. l33t l34l l35l| . The statistical 
errors given by the experiments are smaller than the symbol 
sizes. Systematic errors are not available. 



line shows a linear fit through the data points, while the 
data points 0, El ESI EH are depicted by full symbols. 

At a first glance the energy dependence looks struc- 
tureless. The data seem to follow a linear dependence 
on the beam rapidity y p without any irregularities. How- 
ever, the general trend of the rapidity widths is also well 
reproduced by Landau's model with an equation of state 
with a fixed speed of sound. Nevertheless, there seem 
to be systematic deviations. At low AGS energies and 
at RHIC, the experimental points are generally under- 
predicted by Eq. J2J), while in the SPS energy regime 
Landau's model overpredicts the widths of the rapidity 
distributions. Exactly these deviations from the simple 
Landau picture do allow to gain information on the equa- 
tion of state of the matter produced in the early stage of 
the reaction. By inverting Eq. J2J we can express the 
speed of sound c 2 in the medium as a function of the 
measured width of the rapidity distribution: 



4 ln( x /SNN/2m p ) 
3 



4 ln( % /sNN/2m p ) 
3 Vl 



(3) 



Let us now investigate the energy dependence of the 
sound velocities extracted from the data. Figure [2] shows 
the speed of sound as a function of beam energy for 
central Pb+Pb (Au+Au) reactions as obtained from the 
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FIG. 2: Speed of sound as a function of beam energy for 
central Pb+Pb (Au+Au) reactions as extracted from the data 
using Eq. ||3J. The statistical errors (not shown) are smaller 
than 3%. 



data using Eq. ©. The sound velocities exhibit a clear 
minimum (usually called the softest point) around a 
beam energy of 30A GeV. A localized softening of the 
equation of state is a long predicted signal for the mixed 
phase at the transition energy from hadronic to partonic 
matter pjfjL Wh . l38| . Therefore, we conclude that the mea- 
sured data on the rapidity widths of negatively charged 
pions are indeed compatible with the assumption of the 
onset of deconfinement at the lower SPS energy range. 
However, presently we can not rule out that also an in- 
creased resonance contribution may be the cause of the 
softening |3Sj . 

In conclusion, we have explored the excitation func- 
tions of the rapidity widths of negatively charged pions 
in Pb+Pb (Au+Au) collisions. 

• The rapidity spectra of pions produced in central 
nucleus-nucleus reactions at all investigated ener- 
gies can be well described by single Gaussians. 

• The energy dependence of the width of the pion 
rapidity distribution follows the prediction of Lan- 
dau's hydrodynamical model if a variation of the 
sound velocity is taken into account. 

• The speed of sound excitation function extracted 
from the data has a pronounced minimum (softest 
point) at Ebeam = 3CL4 GeV. 

• This softest point might be due to the formation 
of a mixed phase indicating the onset of deconfine- 
ment at this energy. 
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Further explorations of this energy domain is needed and 
can be done at the future FAIR facility and by CERN- 
SPS and BNL-RHIC experiments. 
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